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Abstract 



I introduce a straightforward technique for the filtering of extended astronomical 

images into components of different spatial scales. For a positive original image, each 

^ ■ component is positive definite, and the sum of all components equals the original 

^3^ . image. In this way, the components are each individually suitable for flux measure- 

fv^ I ments and broadband spectra calculations. I present an illustration of this technique 

O ' on the radio galaxy Cygnus A. 

(^ 

P . 1 Introduction 

One common problem in astronomical image processing is the detection of 
K^ , point or small-scale features against a varying background, or more gener- 

H I ally, the separation of information on different spatial scales. I introduce a 

technique for filtering extended astronomical images based on a suggestion in 
Starck, Murtag and Bijauoi (1998). Images can be decomposed into an arbi- 
trary number of components, each representing structures on different spatial 
scales and whose sum equals the original image. Each component is positive 
definite and therefore can be suitable, e.g., for spectral index or optical color 
studies. The method is similar in some ways to the multi-scale wavelet analysis 
discussed by (Slezak, Durrett & Gerbal 1994). 

I have applied the multiresolution filtering to a number of radio galaxies. I 
find that the spectra of fine-scale features do not follow a simple pattern, 
as opposed to the slow gradients which often characterize the more diffuse 
emission. In this paper, I will show some results on Cygnus A as an example. 
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2 The Method 



For the simplest two-component image separation, the initial image is decom- 
posed into Filtered and Open maps, where Original = Filtered + Open. 

The Open map is made by first replacing each pixel on the original map by the 
minimum value in a box of fixed size which slides across the map, (a process 
called "erosion"), followed by replacing each pixel on the "eroded" map by 
the maximum value in a sliding box across the map. The result is the Open 
map. Then, Filtered = Original — Open. 

The Filtered map now contains only structures with sizes smaller than or on 
the order of the sliding box size. Everything else, by construction, is in the 
Open map. The method can be generalized by decomposing the Open map 
itself into an additional Filtered' and Open' map by operating with a larger 
box on the Open map. This method can be generalized to an arbitrary number 
of steps. In practice, I have found that a change in box size of approximately 
four between steps is necessary to isolate independent components. 

This technique has been implemented in AIPS with a trivial modification to 
the task MWFLT and a RUN file to automate the various steps. It would also 
be quite straightforward to implement, e.g., in IDL. Readers are encouraged 
to contact the author for details. 



3 The Use of Multiresolution Filtering for Spectra 



The measurement of spectra in extragalactic radio sources is fraught with 
problems such as the confusion between overlapping structures, biases due to 
bowls, and the subjectivity of modelling or other techniques used to isolate 
features. Each spectral measurement technique, e.g., direct division, the re- 
gression, or "T-T" method, or even the tomography we have introduced (?), 
has its own particular problems. Multiresolution filtering offers a complemen- 
tary method which can be useful when applied with care. I have done a great 
deal of testing of the robustness of this technique, as will be described in an 
upcoming paper. An illustration of our results on Cygnus A is shown in Figure 
1. 

To calculate the spectral index, I first decompose the maps into two or more 
components. Then, the spectral index can be calculated separately for each 
component, either by straight division or by regression. In cases where there 
is a clear bimodal (or more) distribution of spatial scales, this is a meaningful 
process. If, however, there is a continuous range of spatial scales, then the 



decomposition and the resultant spectra are somewhat arbitrary. They still 
may be useful, however, for seeing whether there are differences between fine- 
and large- scale spectra. 

In general, I find that fine-scale structure has a slightly flatter spectral index 
than the local large-scale structure. For Cygnus A, however, there are locations 
in the source where the opposite is true. The one-dimensional plots at the 
bottom of Figure 1 do not show the full range of spectra, but still give an 
indication of the richness of the fine-scale spectral structure. 

Looking at the spectra of the large scale emission (which might be most appro- 
priate for ageing analyses if we could eliminate the problems due to magnetic 
field variations), we can use this filtering technique to assess/eliminate prob- 
lems with confusion from small scale structures. In the case of Cygnus A, this 
turns out not to be a significant problem, as seen by the close correspondence 
between the spectra of the original (not decomposed) and open maps. 

How can the large-scale and fine-scale spectra be different? First, as seen in 
the numerical simulations such as those presented here by Jones, Tregillis and 
Ryu, there can be filamentary or other structures with very different particle 
acceleration and/or loss histories than neighboring features. Alternatively, the 
local magnetic fields may vary between structures, causing a change in the ob- 
served spectral index even if the relativistic particle populations are identical. 
It is this confusion which causes major problems for standard ageing analyses. 

I have also worked with colleagues to experiment with galactic plane contin- 
uum data and deep optical/IR images. It appears that multiresolution filtering 
will be useful for a wide range of astronomical images. 

This work is supported by NSF grant AST 96-16964 and AST AST 96-19438 
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Fig. 1. Filtering results on Cygnus A. (All based on data from Carilli et al., 1991.) 
The top image is the total intensity at A20cm with a resolution of 1.5". The middle 
image is the open map and the bottom the filtered map with a box size of 4.8" . The 
plots at the bottom show the total intensity in transverse strips along the major 
axis along with three different measures of the spectral index between X20cm and 
X6cm. The x's show the indices from the original maps at two freqeuncies; the open 
circles from the corresponding open maps, and the filled circles from the filtered 
maps. 



